Experimental Test of a New Equality: 
Measuring Heat Dissipation in an Optically Driven Colloidal System 
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Measurement of energy dissipation in small nonequilibrium systems is generally a difficult task. 
Recently, Harada and Sasa [Phys. Rev. Lett. 95, 130602(2005)] derived an equality relating the en- 
ergy dissipation rate to experimentally accessible quantities in nonequilibrium steady states de- 
scribed by the Langevin equation. Here, we show the first experimental test of this new relation in 
an optically driven colloidal system. We find that this equality is validated to a fairly good extent, 
thus the irreversible work of a small system is estimated from readily obtainable quantities. 
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Experimental techniques of the manipulation and ob- 
servation in small systems has developed significantly 
in these years. Recently, new relations in statistical 
mechanics such as the fluctuation theorems (T), d, & 0] 
and the Jarzynski equality]!, 0] have been discovered 
for small systems driven out of equilibrium. These new 
ideas have been experimentally tested @, H, [HI Hlj 
and utilized [12, EH 13 mainly in the systems of macro- 
molecules and colloids using the optical tweezer. 

Recently, Harada and Sasa derived an equality [l5lll6l. 
Il7j relating the rate of energy dissipation to experimen- 
tally accessible quantities described by Langevin equa- 
tions. Although the energy dissipation rate is the fun- 
damental quantity to characterize the nonequilibrium 
steady states, it is generally hard to measure it directly 
in small systems, especially in microrheological systems 
and biological systems such as molecular machines. The 
important aspect of this relation is that it enables us to 
evaluate the energy dissipation rate from experimentally 
accessible quantities such as correlation functions and re- 
sponse functions of fluctuating observables . Moreover, 
this equality was generalized to systems with a memory 
kernelfT8j. In those systems, the fluctuating force has 
a coloured spectrum, not a white one. Such a general- 
ization provides us applications of the equality to much 
wider systems, including rheological systems and biolog- 
ical systems. In this paper, we showed the first exper- 
imental test of the prediction by Harada and Sasa in a 
colloidal system driven by an optical tweezer in a viscous 
fluid. Also, we evaluated the energy dissipation rate as a 
function of frequency using this equality. 

Firstly, we introduce this equality. In the over damped 
case, the Langevin equation is written as 

where x(t) denotes a certain degree of freedom such as 
the position of a Brownian particle. The force F(x(t),t) 



acting on the particle at x(t), can be periodically or ran- 
domly changing in time, or time independent in general. 
sf p (t) is a small probe force to measure the response 
and £(t) is a white Gaussian random force : (£(£)) = 0, 
= 2jk B TS(t — s). In the systems described by 
this Langevin equation, Harada and Sasa derived the fol- 
lowing equality [H, EE QjJ , which relates the ensemble 
average of the energy dissipation rate J(t) to experimen- 
tally accessible quantities as : 



Wo = 7 



C{u) - 2k B TR'(uj) 



Z7T 



where Fourier transform of an arbitrary function A(t) is 
defined as A(u) = A(t) exp(iu;t)dt. A'(lj) denotes 

the real part of A(ui). C{t) is the auto correlation func- 
tion of the velocity : C(t) = (x(t)x(0)) o . R(t) is the lin- 
ear response function of the velocity to a small external 
probe force sf p (t) : 



R(t~ S )fP(s)ds + o(s 2 ), (2) 
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where (• • ■ ) e denotes an ensemble average under the 
probe force f p (t) of order e. v s is the steady state ve- 
locity. Measurement of C(ui) and R'(oj) is usually more 
accessible by experiments than that of J(t). The right 
hand side of Eq. ([T]) vanishes near the equilibrium state 
due to the fluctuation dissipation relation of the first 
kind[H; Q(w) = 2k B TR'(uj). While it has a finite value 
in nonequilibrium states generally. Therefore the integral 
in Eq. (fTJ) corresponds to the extent of the FDR violation. 

We tested this equality (Eq. (JTJ) in an optically driven 
system, in which one small particle was trapped by an 
optical tweezer. The optical tweezer enables us to cre- 
ate a potential for small particles [20j. Nonequilibrium 
steady states(NESS) can be established by introducing 
persistent forcing to the Brownian particle. 

The experimental setup is shown in Fig.[T] Infrared 
laser(Spectra Physics, 1064nm) was introduced to the mi- 
croscope(01ympus, 1X71) and focused by an objective 
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lens(01ympus UPLSAPO,100x NA=1.4) to create an op- 
tical tweezer. One carboxylated polystyrene particle of 
a diameter of 0.984±0.023^m(Polysciences,#08226) was 
trapped in distilled water. Trajectory of the particle was 
captured by a CCD camera(Hamamatsu, HISCA C6770) 
with 472 [frames/s] and analyzed to know the particle po- 
sition. The temperature was kept at 27.0±0.1[°C] by 
water circulation around the objective lens. 
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FIG. 1: Experimental Setup 

To achieve nonequilibrium steady states(NESS), we 
switched the position of the tweezer between two sites 
±L/2(L = 460.8 ± 0.14[nm]) at random according to 
the Poissonian process by tilting the angle of a Piezo 
mounted mirror(PI, S-226) controlled by a PC. The cor- 
responding Langevin equation becomes 



7±(i) = — k 



x(t) - -&(t) 



where k is the stiffness of the tweezer. x(t) is the hori- 
zontal coordinate of the particle position (see Fig.[^a)). 
a(t) = { — 1,1} is a random variable obeying Poissonian 
process with a transition rate A; (a(t)a(s)) = cxp(— 2A|t— 



,s|). The particle is sufficiently small (m/7 
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thus the over damped assumption is supposed to be valid. 

At first, we trapped a particle with fixing the po- 
sition of the tweezer and obtained a histogram of the 
particle location(l,703,884 points). Then, according to 
the Boltzmann distribution, we calculated the potential 
energy of the particle, which was well fitted by a har- 
monic potential with a spring constant of k — 117.5 ± 
0.48[k B T//im 2 ] = 0.4870 ± 0.0020 [pN//im]. 

To know the response function we shifted the 

position of the tweezer at t = with a distance of 
d =190.2±0.11[nm] (Fig.[3|) and tracked particle trajec- 
tories in the relaxation processes in the equilibrium 
state(EQ) and four NESSs(the mean switching rates were 
varied from A=1.0 to 4.0 with a step of 1.0[Hz]). The 
corresponding probe force ef p (t) becomes kd(t > 0), 
(t < 0). Figured] shows that the ensemble aver- 
aged trajectories of the relaxation processes in EQ and 
NESSs coincided well. We fitted an exponential func- 
tion {x(t)) jd = 1 — exp(— Lu c t) to the trajectory in 
EQ. The cutoff frequency was estimated as w c =27.61± 
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FIG. 2: Switching: (a)To achieve nonequilibrium steady 
states(NESS), position of the optical tweezer was switched be- 
tween two sites(±L/2) using a Piezo mounted mirror at ran- 
dom with Poissonian process. In the equilibrium state(EQ), 
the position of the tweezer was fixed. In NESS, we decided 
whether to switch or not at random every 4ms. One switch- 
ing takes 5[ms] due to the mechanical response of the mirror 
unit. (b)Typical trajectory for EQ and NESS(A = 1.0). Black 
line is the trajectory of the optical tweezer. Gray line is the 
trajectory of the particle. L — 460.8 ± 0.14[nm]. 
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FIG. 3: Measurement of Response Function : (a)We 
shifted the tweezer position with a distance of d = 
190.2±0.11[nm] at t = and observed the relaxation processes 
to the new laser position. We repeated this and obtained an 
ensemble average. (b)Similarly in NESSs, we added a bias d 
and observed the relaxation processes. 



0.05[rad/sec], which yields the friction coefficient 7 = 
k/uj c = 1.763 x 10- 8 [kg/s][2l|. 

From the fitting curve in Fig. 2] we obtained the re- 
sponse function R'(li) = uj 2 /(u 2 + uj 2 )^. We compared 
this response function R' (uj) and correlation functions of 
the velocity C(ui) obtained from the trajectories of the 
particle. In Fig. [5] we found a good coincidence between 
C(uj) and 2k-QTR'(uj) in the equilibrium state. Thus, the 
fluctuation dissipation relation was confirmed. While, 
these had differences in nonequilibrium steady states. 
Also, the higher the switching rate A was, the larger the 
deviation of C(ui) from 2ksT R' (w) was. 

Equation |l| implies that the product of 7 and the 
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FIG. 4: (Color) Response Function : The trajectories of 
the relaxation processes in EQ (equilibrium state(no switch- 
ing)) and NESSs(nonequilibrium steady states) with mean 
switching rates A= 1.0, 2.0, 3.0, and 4.0[Hz]. Dashed 
line is the fitting curve : {x(t)) /d = 1 — exp(— uj c t) for 
EQ. d was 190.2±0.11[nm] and the cutoff frequency lu c was 
27.61±0.05[rad/sec]. We averaged 7,854 trajectories for EQ, 
A = 1.0,3.0, and 4.0 respectively and 8,085 trajectories for 
A = 2.0. 
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FIG. 5: (Color) C{lo) and 2k B TR'{uj) : Dashed line is 
the response function 2k B TR'(uj) = 2k B Tu 2 /(uj 2 + lo 2 )j. 
52[runs]x69.17[s] (l,703,884[points]) in EQ and NESS at A = 
3.0, 4.0 respectively and 53[runs] x69.17[s] (l,736,651[points]) 
in NESS at A = 1.0 and 2.0 respectively were averaged for 
C(uo). Moving averages were taken in to direction. 
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area between C{uf) and 2k-QTR'(Lo) in Fig. [5] corresponds 
to the energy dissipation rate. The irreversible energy 
dissipation should be balanced with the energy input in 
steady states, thus (J) = (F(x(t),t) o y(t)), where o de- 
notes the Stratonovich multiplication 23] . On the con- 
trary, in our setup it is possible to evaluate this (J) be- 
cause the force exerted by the optical tweezer can be cal- 
culated from the displacement between the particle and 
the center of the tweezer as F(x(t), t) — —k{x{t) ~ -§c(i)) 
(Note that this is not always possible in other systems.). 
We compared these two values in Fig.[S]Ja) and confirmed 
that the equality is valid to good extent within statisti- 
cal error. In one switching, the particle is supposed to 



FIG. 6: (Color) (a)The left hand side(LHS) and the right 
hand side(RHS) of the equality{T} were plotted against the 
mean switching rate A [Hz]. For RHS, we integrated only in 
the region w < 850 because the difference between LHS and 
RHS should vanish at sufficiently high frequencies. Dashed 
line is the theoretical curve [22]. The standard errors of RHS 
were around 1[]£bT/s] at all switching rates, (b) Validation of 
the equality for each w(Eq. [31; 7(w) (blue line) vs 7[C l (w) — 
2&bT-R' (w)] (black line). We used the same data for C(w) 
and R'{u) in Fig.[S]to calculate j[C(u) - 2k B TR'(uj)] 
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travel a distance of L on an average. This corresponds 
to an energy injection of kL 2 /2 = 12.5 [kaT] per one 
switching and XkL 2 /2 = 12.5A[keT] per unit time. How- 
ever, insufficient relaxation to the switched position due 
to fast switchings causes a deviation from this value. The 
dashed line in Fig.[6ja) is the theoretical value including 
this effect [22 1 . 

It is possible to derive a local equality (Eq. (TT|)): 



I(lu) = 7 C{uj) - 2k B TR'(uj) 



(3) 



where I(u>) is the Fourier transform of the cross correla- 
tion of the force F(x(t),t) and velocity v(t) at different 
times: 

I(t) ee [(F(x(t),t) o v(Q)) + (F(x(0),0) o v(t)) ] 12. 

Equation (J3]) implies that the equality holds at each 
a; [la]. This local relation provides us a chance for the 
more detailed experimental test of the equality. Fig- 
ure EJb) shows a good agreement of the both sides of 
Eq. ([3]) at each to. As we increased the mean switching 
rate A, the frequency and height of the peak of the local 
dissipation rate increased. 

We showed the first experimental demonstration of the 
Harada and Sasa's equality. Using this equality, we eval- 
uated the rate of the energy dissipation for an optically 



driven particle with a precision of around lfksT/s] from 
the correlation function and response function of the par- 
ticle velocity. Note that this equality is not applicable 
to systems with retarded friction, since it is based on 
a Markovnian Langevein equation. However, systems 
which are more interesting and practical often have a 
retarded friction. For example, molecular machines are 
working in extremely crowded environment in cells. In 
such a situation, non-Markovnian behaviour is expected. 
Also, microrheology is attracting extensive studies re- 
cently. As we described, Deutsch and Narayan recently 
generalized the equality to these situations, namely sys- 
tems described by a generalized Langevin equation [18[. 
If this generalization is confirmed by experiments, it pro- 
vides us applications of this equality to much wider sys- 
tems. This will be the primary focus of the forthcoming 
experiments. 
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